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NATIONAL ADVISORY COMMITTEE FOR AERONAUTICS 
TECHNICAL MEMORANDUM 1264 


CALCULATION OF THE LATERAL-DYNAMIC STABILITY OF AIRCRAFT* 
By А. Raikh 


The object of this report is to present a method of computing 
the lateral-dynamic stability of an airplane. 


Graphs and formulas are given with the aid of which all the 
aerodynamic coefficients required for computing the lateral dynamic 
stability can be determined, А number of numerical examples are 
given for obtaining the stability derivatives and solving the 
characteristic-stability equation. Approximate formulas are derived 
with the aid of which rapid preliminary computations may be made and 
the stability coefficients corrected for certain modifications of 
the airplane, А derivation of the lateral-dynamic-stability equa- 
tions is included. 


| 


— 


INTRODUCTION жеп 


In the present stage of development ої airplane design, con- 
stantly increasing requirements are imposed on the airplane as 
regards stability and maneuverability. if the area of the vertical 
tail surface and the transverse dihedral of the wing are unfavorably 
chosen, the airplane will be subjected to large lateral motions at 


the least gust of wind and the pilot will be under the constant --- 


necessity ої applying the controls. Some airplanes because of their 

flying qualities are generally approved by pilots, Computations show 
that such airplanes are spirally unstable only at sufficiently large 

angles of attack but are spirally stable at small angles of attack, 


A method of computation of the lateral-dynamic atability is 
presented herein. The computation is illustrated by examples of 
the computation of a Northrup 2E airplane, which has received the 
highest commendation of pilots for its piloting characteristics with 
regard to stability and maneuverability. 

The difficulty of computing the dynamic stability lies in the 
fact that a large number of aerodynamic coefficients, the so-called 
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rotary derivatives, are required for the computation. These deriva- 
tives are not encountered in the usual aerodynamic computation. In 
order to facilitate the computation, a number of curves are given 
with the aid of which these magnitudes are rapidly and simply deter- 
mined, The numerical examples given for computing the coefficients 
further simplify the problem, 


A comparison of the aerodynamio coefficient obtained by compu- 
tation from graphs and that determined from tests in the wind tunnel 
is presented herein, In order to check the accuracy of the final 
data, a comparison is given of the period of oscillation ani the 
time taken for the damping of the lateral motion as obtained from 
flight tests and by computation, respectively. The comparison shows 
that the computation gives good results. 


Notwithstanding the apparent complexity, the entire computation 
on the lateral stability from the available aerodynamic coefficients 
takes no longer than 1.5 to 2 hours. The author expresses his thanks 
to А, I. Silman, a candidate for a technical degree, for a number of 
valuable suggestions utilized in this work, 


SYSTEM OF AXES AND NOTATION 


The system of axes is shown in figure 1. The origin of the 
coordinates is at the center of gravity ої the airplane, 


The X-axis ig in the plane of symmetry of the airplane and dir- 
ected parallel to the velocity in steady flight. The Y-axis is рег- 
pendicular to the plane of symmetry of the airplane and directed to 
the right of the pilot. The Z-axis ів in the plane of symmetry рег- 
pendicular to the X-axis and directed downwards. | 


тҮ? forces along corresponding axes; positive if 
direction agrees with direction of axes 


L, M, N moments about corresponding axes: moment L is 
positive if it causes right wing to drop; 
moment M is positive if it causes tail to 
drop; moment N is positive if it causes 
backward motion of right wing | 


u, Y, W projections of linear velocity on corresponding 
axes; positive if direction of projections 
agrees with directions of axes а 

р, 4, г projections of angular velocities; have ваше posi- 


tive directions as moments L, M, N 


-— = 
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Ф, 9, У 


А, В, С 


angles made by given direction with Х-, Ү-, 2- 
axes respectively; have same positive direc- 
tions as L, М, N 


moments of inertia about Х-, Ү-, Z-axes 
the centrifugal moment of inertia (Еа. note, 
product of inertia with respect to Х- ага Y- 


axis 


angle of attack, degrees 


angle of attack computed from line of zero lift, 


degrees 


sideslip angle, В = -arc віп v/V, radians (МАСА 
comment: 8 as used in this report is of 
opposite sign to В as used in American and 
British reports,| 


angle of inclination of Plight path to horizon- 
tal (positive for lift), degrees 


angle of dihedral of wing, angle between plane 
of chords and plane at right angles to plans 
of symmetry and passing through chord at Бір, 
degrees 


angle of sweepback, angle between focal line at 
0.25 chord computed from leading edge and 
plane perpendicular to axis of fuselage, 
degrees 


angle between line joining center of gravity 
of airplane with geometric center of verti- 
са] tail surface (with part of fuselage) 
and line of zero lift, degrees 

area of wing 


area of vertical tail with part об fuselage 
(fig. 11) 


area, of lateral projection of fuselage 


span 








4 
29 
^ 
À = 52/8 
1 = og/ce 
inf 
it 
ір 
б 
E 
Y 
NER 
y 
ота 
2 
бі = 
! wu 
2 
N 
С. - 
п 2 
ev 
2- Sb 
t 
Or =. z 








МАСА TM 1264 


chord at root of wing 

chord at tip of wing, effective tip chord is 
determined by prolonging lines of leading 
and. trailing edges 

aspect ratio 


taper of wing 


distance from nose of fuselage to center of 
gravity of airplane 


distance from rudder hinge to center of gravity 
length of fuselage . 
density of the air 


coefficient of effectiveness of vertical tail 
surface 


coefficient of lateral force 


coefficient of rolling momenti  . 


coefficient of yawing moment 
117% force coefficlent 


nomimensional derivative of lateral force with 
respect to sideslip angle (8 in radians) 


nondimensional derivative of rolling moment with 
respect to sideslip angle 
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3c | 
ng See nonüimensional derivative of yawing moment with — 
98 | respect to sideslip angle 
ӘС д в о 
ір = — noniimensional derivative of rolling moment with > U 
а 25 respect to angular velocity of roll ZEE 
дс 
Ny ша nondimensional derivative of yawing moment with 
ә Бе respect to angular velocity of roll | К 
ly өрі noniimensional derivative of rolling moment with 
ә 55 respect to angular velocity of yaw 8 fons 
Dy = 901 nondimensional derivative ОҒ yawing moment with 
3 P respect to angular velocity of yaw 
ід = -$ coefficient of moment of inertia about X-axis 
mb 
іс - 25 coefficient of moment of inertia about Z-axis 
mb | 
а weight of airplane i 


masB of airplane 


relative density of airplane | AE т. 


€ 
Є 
.2m 
oSb 

T = 27 unit of time 
z nondimensional time | ал 
Аз, 


Ao, Ау, Ад coefficients of characteristic equation 
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і, EQUATIONS OF LATERAL-DYNAMIC STABILITY OF AIRPLANE 


If the rectilinear steady flight of an airplane is considered , 
the motion of airplane after a certain disturbance is described by 
a system of six equations with віх variables: 


u, Wy 4; and B, р, r 
If. only small deviations of the airplane from steady recti- 
linear flight are considered, it may be shown (reference 1) that 
the system of віх differential equations breaks up into two ідйе- 


pendent systems of equations, In the first system of differential 
equations the following variables enter: 


u projection of linear velocity on X-axis 
w projection of linear velocity on Z-axis 
а projection of angular velocity on Y-axis 
То the first system of equations, there enter only variables 


characterizing the motion of the airplane in the plane of symmetry 
(that is in the plane Х02). In the second group of equations there 


enter the variables: 


В angle of sideslip, determined by projection of velocity оп 
Y-axis (sideslip velocity т) 


p projection of angular velocity on X-axis. 
r projection of angular velocity on Z-axis 


The variables В, p, and r characterize the lateral motion of the 
airplane, 


By assuming the nondependence of the lateral motion of the air- 
plane on its longitudinal motion, the equations of the lateral дупа- 
mic stability can be derived, 


Yor steady voctilinsar flight without sideslip 
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Х-Ү =Z =L «М = ХК = 0 


(1) 


М = У 
where У is the velocity о? the airplane. 
After a certain lateral disturbance, the magnitudes р, r, and v 


will vary with times, whereas the magnitudes q, м, and u may be con- 
sidered constant, | 


Because at the initial instant р - Г - 7 - 0, the disturbances 
of the magnitudes р, г, and у are ИН by the same symbols: 
Ар = р - ee. te 
Ат = r MEE 
АУ = У 


Evidently, the lateral motion of the airplane will be characterized 
by three equations: (a) the equation of the equilibrium of the for- 
ces along the Y-axis, (b) the equation of the equilibrium of the | ТЩ 
moments of the forces about the X-axis, and (с) the equation of the і 
equilibrium of the forces about the Z-axis, Ву the fundamental 
equation of mechanics, the product of the mass of the airplane by 
the projection of the absolute acceleration on the Y-axis is equal 
to the projection of all external forces acting on the n on 
the Y-axis, | 


In the general case, the projection ої the —ÓÀ accelera- 
tion of the Y-axis is equal to (reference 2) 


dv | TP 
ЗЕ t Tú ~ РУ | | n 


and the equation of equilibrium of the forces адна the Y-axis will | 
be of the form на бац 4, ы 


а (@ +0) = | (2) 
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For the case under considerationi 
u= У етй Ww = О ' aaa 


In equation (2), У denotes the sum of the projections on the Y-axis 
of all the external forces acting. on the airplane so that 


Y = Y, + tw (5) 


where 
Ха projection of aerodynamic force on Y-axis 


У, projection of weight on Y-axis 


In what follows, it is assumed that the aerodynamic force acting 
on the airplane during a disturbed motion and the small deviations 
from the condition of steady rectilinear flight depend linearly on 
the disturbances, and because at the initial instant no lateral force 
acted on the airplane the aerodynamic force is given by 


Ха . Sr у + s: р + xr 
ТР the small magnitudes EH p and S r are neglected, the fol- 
lowing equation is obtained: | 


Ха = S v А (4) 


In symmetrical flight there аге no gravitational forces acting on 
the airplane along the Y-axis, If the airplane rolls, a component 
appears on the Y-axis. When small oscillations occur, the force 
acting along the Y-axis for a rolling angle ф is equal to 

(figs. 2 and 3) 


— 


G сов 9“ф 
—raasaan ИИИ л rr rc 


lin а rigorous derivation, the term pw drops eu ав а mag- 
nitude ої second order of smallness. 
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where 0 із the angle of inclination of the flight path, Similarly, 
for the rotation of the airplane about the Z-axis by an angle Y 
the component along the Y-axis (fig. 4) is 

С sin 6y 
Hence, for small deviations 


Y, = G cos 0$ + G sin 6-ұ = @ cos 6(Ф + ten Oy) 


But 
С сов 8 = ст, EY. 
Therefore, 
уа 
Үр = Cr, B2 8(ф + tan 6. у) | (5) 


For small oscillations ої the airplane, it may be assumed that 
У = - V sin pw - УВ 


When this relation is taken into account and equations (3), {аў, 
and (5) are substituted in equation (2) 


о 
a (- 8 + =) вас Fs + tan өзү) (6) 


The equations of equilibrium of the moments about the Х- and 7-ахев 
are of the form (reference 3) 


Ар ра. 
З at " dt L те 
4. й | i I 
= - с.а | 
С Е at N E Ж 


The moments of the aerodynamic forces Т and М оп the free motion? 
of the airplane after a disturbance will depend on the angular veloc- 
ities р and г and the angle of sideslip В, For the small dis- 
turbances р, г, and В, the aerodynamic moments тау be represented 
as follows by a linear dependence on р, г, and В. 


—— — s TI 
ZThe free motion of the airplane after а disturbance with locked 
controls, 
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à 
р. р. Ў: 


Ў Шы Шы 


Equations (7) now become 


a ат а, aL ӘТ, 
AH - E S£ = Š + РЭ: 
(8) 
dr ƏN ON 


dp -ƏN 
ба - Ча = 36 Ë + S P + г 
ТР equation (6) is differentiated together with equations (8), а 


system of three linear homogeneous differential equations isa 
obtained because p = 90/4% and r = dy/dt. 


2 2 
av (- 8.22) - 3118 or ӨЙ slp + = tan 6) 


ate at 
а ду дт, ОТ, aL 

А XE - Е +£ = Эв B + 55 р + 5 T (9) 
dr d ƏN ƏN ƏN 

C= - ESS за + S; P + S " 


Nondimensional magnitudes are now introduced by the equations 


2 
ev 


! (10) 
oV 
N = Cn 727 Sb 
mbé 
dc Aou 





Р | | | 
Су is here the lateral force coefficient which should not be 
confused with the previously used notation for the Lift coefficient. 
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mb? 
indo r 
2 
E = ip” 
= рь 
P = зу 
z. rb 
“аў 
ac 2 2 ac 
S = З-за Jp = Зв 
ac 
ac 2 2 әс 
ОН. сд SPS so = ng SE sb пр = xs 
aL у, ры - 90; сү? pve 3c 
3p Р = уру" a SP = Ply 2 Sb ір“ ЗЬ 
2% 27 27 
2 әс 
SE к = Flp E Sb ly = — 
9 37 
" 2 
My ба, SP в is: 
2v 
д ~ pve On 
ЗЕЕ = Го» Ё5- Sb “е = ть 
2v 


ТР all terms on the left in equation (9) are transferred and 
the nondimensional magnitudes are introduced, the system becomes 


12 
42 ат dr gv? 4 оў? 
шу av b att 2. 2 Syg sË 0, S "n S! 
mb? 2V dp , mb“? гу dr ру? 
AT dt "E 4 45 3 
4 b at 4 bat” 2 
2 Е 2 С 
, me 2v dr | m 2v ар ру 
C 4 96-18 4 b dt- 8 


ау ы 
= (p+r tan Ө) 


2 
pv" gr 
Sb Фр == вых - | SbIgR. 
2 2 
- pv _ ov 
Зра Р - 75 Sbn,r- = SbngB 
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= 0 
NI 
C 
2:0 M 
= Q 


If the first equation of the system is divided Бу о2У582/п, 


2 
the second by 1, 25°, 





equations are obtained: 


m \ ав 1 ap m r 2m- 
pSV} 4.2 * 2 УВ dt ЗҮ + 2 ЫР” 
CL 
2 собр tan 9Хг = 0 
ЇВ ошар lp- E m аг 
~ 14 P + coy dt - i4 P “ i, 057 dt 
n _ 2р- 18 m dp m dr 
“3g P - Tç P - Я, 057 dt + 057 dt 
By setting 
ты 
psv 
_2m 
В = 57 


and the third by i, 


" 
зіс 


pV^Sb 
2 





, the following 


(11) 


(12) 
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Moreover, if the signs in the second and third equations are changed, 
the following system of equations in nondimensional form is obtained: 


` 2 C w € 
> ав 1. ав, CL dr | CL = = 
3 ag гаў НРУ Bur LP BSc 
i = 1 — 1 | 
в ёр. ps Ear м. (аз) 
Т B Fa кү + 1 г 0 
А а% А А 4% А 


30, дс, дс, 
The magnitudes Jg => ng =p and їв = 35 are determined 


from the usual tests in the wind tunnel and are denoted as the sta- 
tic derivatives, The magnitudes ір, lrs Пр; and ny characterize 


the moments of the aerodynamic forces in rotation and for this rea-. 
gon they are known as the rotary derivatives, 


2, THEORETICAL DETERMINATION OF DERIVATIVES OF 
LATERAL DYNAMIC STABILITY 


As has been pointed out, all aerodynamic coefficients are 
referred to a system of axes: the X-axis is parallel to the veloc- 
ity and passes through the center of gravity of the airplane at the 
initial instant; the Y-axis passes along the span of the wing to 
the right of the pilot; and the Z-axis is perpendicular to the first 
two and directed downward?, The positive direction of the angles of 
rotetion and of the moments are assumed to be the following: 


about the X-axis, motion from the Y-axis to the Z-axis a 
about the Y-axis, motion from the Z-axis to the X-axis "р un. 
about the Z-axis, motion from the X-axis to the Y-axis | | 


This system of axes is fixed in the airplane so that the X-axis 
coincides with the velocity of the airplanes only at the initial 





“This system of axes is in accordance with the terminology 
adopted by САНІ, May 1959, 
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instant of time. Рог a disturbed motion of the airplane, the X-axis 
of the fixed system of axes makes an angle В, which is termed the 
sideslip angle, with the flight velocity. 


The moments about the X- and Z-axis are denoted by L and N, 
respectively, and the angular velocities about the X- and Z-axis Uy 
p and r, respectively. 


ROTARY DERIVATIVE by 


In the following discussion, all the derivatives in the below- 
stall range will be considered. Іп the case of rotation of the air- 
plane about the X-axis (angular velocity р), an aerodynamic damping 
moment arises about the same axis, This moment may be determined by 
the equation : 


220) 97 5 


“аў p 2 
where 
Ly moment coefficient 
P> = nondimensional velocity of roll. - 


Неге as usual (reference 1), а linear dependence of the aerodynamic 
reaction on the disturbance is assumed, The rolling moment is pro- 
portional te the fourth power of the linear dimension, and therefore 
the moment on the tail surface for an angular velocity p is neg- 
ligibly small as compared with the moment acting on the wing. 

Hence, in computing Ly only the moment acting on the wing shall 


be taken into consideration, 


On rotation of the airplane about the longitudinal X-axis with 
angular velocity р? 0 (the right wing drops); at a flight veloc- . 
ity V, it ів found that the angle of attack Increases on the right 
wing. А% the distance y from the plane of symmetry 
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Similarly, on the left side of the wing the angle of attack decreases. 
A moment therefore appears opposing the rotation. Thus, for р» 0, 
L « 0; and the coefficient of the derivative ої the rolling moment 
with respect to the angular velocity of roll is always (at below- 
stall angles of attack) negative (fig. 5) 1р < 0. Figure 6 gives 
curves by which 12 for wings of various aspect ratios and tapers 


may be determined (references 4 and 5). 
Example: п = 1.94, 6734. ЕИ 


From figure 6: 1. = -0.465. 


р 
it сап be seen from the curves that by increases in absolute 
value with an increasing aspect ratio of the wings and that ly 
decreases with increasing taper 1, | | 
The value ої by is computed on the basis of the general theory 


of the wing with account taken of the change in the distribution of 

the circulation along the span due to rotation without considering 

the change in velocity along the span. At below-stalling angles, : 
the derivative by theoretically does not depend on the initial 


twist of the wing nor on the angle of attack of the entire wing. 
The magnitude 1р varies little for different wings, as шау be seen 
in figure 6, The mean value is ір = -0.45, The sign of by may | 


change only in the above-stalling range in the region of autorotation 
of the wing. s de 


ROTARY DERIVATIVE n 


In the rotation of the airplane about the X-axis, the angle of 
attack varies along the wing span and аз а result the lift force 
and the induced drag will not be the same at symmetrical elements. 
Hence, when the airplane rotates about the X-axis, the airplane will 
be acted on by a moment of the aerodynamic forces both about the 
X-axis (15) and &bout the Z-axis (пр). Ав before, consider that 


the moment &bout the Z-axis arising for angular velocity p is pro- 
portional to the angular velocity of rotation, so that 
_ ро 
N = ру MpaSb 


where q is the dynamic pressure, 
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The moment of yaw due to the angular velocity depends Little on 
the tail surface so that in the computations for the dynamic stabil- 
ity only the moment due to the wing is taken into account. 


For positive rotation in roll p> 0 (right side of wing drops), 
the results show that on elements of the right and left wing the lift 
force vector L and the vector of induced drag D are rotated and 
changed (fig. 7}. The airplane will therefore be acted on by a 
moment about the Z-axis (moment of yaw), which for negative angles 
may be positive and with increase in the angle of attack of the wing 
varies linearly so that at moderate angles of attack np is always 


negative (ny < о (Ріє. 7)). Curves of the values пор Op (refer- 


епсе 4) for untwisted wings* of various tapers and aspect ratios are 
shown in figure 8, The curves give the values of Np/ da where о, 


із the angle of attack in degrees and is computed from the line of 
zero lift, Hence, in order to actually compute the value of np, it 


is necessary to multiply the value obtained from the curves for the 
given wing (given aspect ratio and taper) by the angle of attack (in 
degrees) computed from the line of zero lift. . 


For example,.n = 1.94, A = 6,54, а, = 99 (а, is computed from 
the angle of zero lift). ре figure 8, nj/a, = 70. 00595 алд, there- 
fore, (пра = 9? = -0.00395X9 = -0.0556. 


ROTARY DERIVATIVE 1, 


For positive angular velocity r, the left wing moves forward, 
that is, the left wing has a greater velocity than the right wing. 
Therefore, the lift force on an element of the left wing will be 
greater than the lift force on the symmetrical element of the right 
wing and as a result the airplane will be acted on by an aerodynamic 
moment tending to lift the left wing, that is, a positive moment 
about the X-axis. With increase in the angular velocity r, the 
moment L increases in magnitude and the derivative of this moment 
with respect to the angular velocity Т = rb/2V will be positive 
(at positive 1)“. 





“Тре value of np little affects the characteristic of the 
dynamic stability and therefore & correction for the twist of the 
wing need not be introduced. 


Here, as. previously, linear dependence of the moment on the 
angular velocity ів assumed as is permissible because only small 
angular velocities are considered and the terms containing the 
squares of these velocities may be neglected. | 
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rb 
Le 27 1,а8р 


The coefficient 2, is determined Гог the wing and is propor- 


tional to the angle of attack аа, Evidently, both the wing taper 
and the aspect ratio will affect the value of 1... Curves of the 


variation of 1,/a, for untwisted wings of various tapers and 


aspect ratios are given in figure 10. From the curves it may be 
seen that 1, decreases with increase in the taper and increases 


with increase in the aspect ratio. 


It is also apparent in figure 10 that 1,/а, for various types 


of airplane changes little with wings of different tapers and aspect 
ratios, The mean value is 2,/a, = 0.017 
0.018 


In order to obtain the value of 1, at а given angle of attack, 
the value Т/с obtained from the curves must be multiplied by the 


angle of attack (in degrees) for which the computation is conducted 
(the angle of attack is measured from the line of zero lift). 


For example, m = 1.94, À = 6.54, a, = 9°. From figure 10 
1-/0, = 0.0185; 1, = 0.0185Х9 = 0.1665. The vertical tail surface, "- 
will also affect the derivative 1, of the airplane. The part of о. н 
due to the tail may be obtained from the equation ----- 


і, 
2 I 
1 5 
„= КД) “ ein 2(х- 8), (14) 
whero 


a, = (9Ст/да), (fig. 11). 


K+ empirical interference coefficient at tail, mean 
value = 0.8 E 22 
т, distance from center ої gravity of airplane to 


ruider hinge | NE CC 
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ВЕР area of vertical tail surface with part of fuselage 
b wing span 
8 wing area 
x angle between line of zero Lift and line connecting 


center of gravity of airplane with geometric сеп- 
ter of vertical tail surface (including part of 
fuselage) (for a low wing airplane x > 0) 


da angle of attack (from line of zero lift), degrees 
‚ s СЁ Ser 
Х = 12° | 


To find 1, for a, = 90, from equation 14у - — 


lng = 2.05 х0.8(0.589)%х 0.093 X 0.105 = 0.00242 


CORRECTION ON DERIVATIVE 12, FOR TWIST OF WING 


In order to introduce a correction for twist, it is necessary 
to use figure 12. Only the practical determination of the correc- 
tion will be considered here The correction for twist for wings 
of various aspect ratios and the taper чүш 3 can be found by use 
of figure 12, ы: 


For а taper different from 3, the correction would be somewhat 
different but the error will be very small and it may therefore be 
assumed that the correction for twist does not depend on the wing 
taper. Assume, for example, the correction for the wing is found 
to be A= 6. 34, the twist of the wing starts at the distance 
0.226 b/2 from the root, increases linearly, and at the tip attains 
the value 1,50, 


At the distance y = 0.226 5/2 from the root, the twist is 
equal to zero (fig, 13(b)). In figure 13(a) at y- 0.226 Ъ/2, 
A, = 0.001. Subtracting this value from D (fig. 13(a)), 0.0168 is 


SThis method and the curves are taken from reference 5, where 
a more detailed description of the curves is found, 
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obtained, This number is plotted in figure 15(с) and corresponds to 


an angle of twist equal to zero; therefore, on the curve ої fig- 
ure 13(c) the point (0, 0.0168) is plotted. It can be seen that at 
the distance 0.4 b/2 the angle of twist is 0.3° (fig. 13(b)) and 
from figure 15(а) Аг = 0.005; l/a = D - Аҙ = 0.0158. In fig- 


ure 13(с) 0.39 is plotted ав the abscissa and -0.0158 as the ordi- 
пабе, and so forth.. A certain area is then obtained on the plani- 
metering for which the correction in 2, is obtained, 

For the given case, ly twist, = 0.0144. 


a ee ee asa not depend on the angle of attack of the 
wing’. 


For the total value of 1. 
ly 
by 42 ба + Түр + br twist (15) 


Example: find 2,: A= 6.34, n = 1.94, а, = 99, а. = 2.05, 
К. = 0.8, 1,/5 = 0.389, 5,./8 = 0.095, x = 129 


The twist starts at 0,226 of the half span, increases linearly, and 
at the tip reaches the value 1.59 


ln = (1,./а,)а, + lpg + ly twist = 0.1665 + 0.00242 + 0.0144 = 0.1833 ` 


ROTARY DERIVATIVE п, 


During an angular velocity of rotation about the Z-axis (angu- 
lar velocity of yaw г) a yawing moment arises, which depends оп 
the vertical tail surface, the fuselage, and the wing, At an angu- 
lar velocity r, on account of the difference in the velocities at 
the symmetrical "elements of the wing, different drags аге obtained 
ав a result of which there arises a moment about the Z-axis (moment 
of yaw N). Evidently, the moment М due to the wing will be neg- 
ative for positive values of т (fig. 14). Similarly, for the ` 
angular velocity г» 0 a negative yawing moment is found to ss 
because of the fuselage and the tail. 





‘The value of 1 һ is positive for the case where the 


r twis 


angle of attack at the tip ої the wing is greater than чне angle of ` 
attack at the root section, t 


. 
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The yawing moment will increase in absolute value with increase А 
in г; hence, 


rb 
N = 57 пуид80 : i 


пу < 0 (fig. 14). 


The effect of the tail and the wing on the derivative п, will now 
be considered, 


Effect of Vertical Tail Surface on ть 


Tn the oscillation of the airplane about the Z-axis, that is, at 
an angular velocity r, the tail surface is found to be acted upon by 
the force 


2 : 
аст, гі ОЗУ 5 
«ПОЛ СП Ў сов(х- в 
? (за), ү 2 в) 


where х ів’ %һе angle between the line of zero lift and the line 
connecting the center of gravity with the geometric center of the 
vertical tail surface with part of the fuselage. For a low-wing 
airplane, the angle x > 0, Іп the expression for the force during _ 
oscillations abowt the Z-axis there enters сов (х-о,) because the 


system of axes is fixed, If the derivative of the moment of this 
force is taken about the center of gravity with respect to the non- 
dimensional angular velocity Y = rvb/2V, 


2 [асу 1% дб 3 
Mee = 2 (да јь Б ш 998 (Х- ча) 


Finally, the nondimensional derivative п, is obtained by dividing "t 
this expression by  oSV^b/2 


ас 1,12 8 у ла 
мя (ве ate 


In this equation no account is taken of the effect ої the fuselage 
on the flow at the vertical tail surface, The indeterminancy of 

the magnitude of the interference at the vertical tail. surface makes 
the theoretical computation of the moments due to the vertical tail 
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surface inaccurate. An empirical coefficient K, characterizing 


the velocity interference at tho tail surface is therefore intro- 
duced. Moreover, setting | 


gives 


By permitting а small error, it may be assumed that cog^(x - сы) 21 


and then 
2 
А Ste 
ыға am sem а. (=) = | (16) 
where 


Uy distance from center of gravity of airplane to rudder hinge 
р wing span | | 
бір area of vertical tail surface with part of fuselage 
5 area of wing 

For the coefficient К, the value 0.8 may be taken. 


The value of ақ depends on the aspect ratio of the vertical 
tail surface. | 


Figure 11 shows the dependence of a, on the experimentally 
Obtained aspect ratio A, (reference 6). The mean value of the 
coefficient а. for airplanes of the usual type is 


a+ = Seo 


In order to compute the aspect-ratio of the vertical tail surface 


À,, the span of the tail b, must be taken as shown in figure 11. 


In the tail surface area Sip there is included the end of the 
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fuselage (hatched area in fig. 11). The effect of the fuselage is Š 
not taken into account in equation (16). Іп the same manner ав Гог 
determining the damping of longitudinal oscillations, the coeffi- 
cient 1.25 in equation (16) ів introduced бо take account of the 
fuselage. Then 


2074 


1\2 8 
t| Mtr 
Depp = а.а | +) “а (16a) 
Example: K+ = 0.8 аъ = 2.05 | 
t,/b = 0.389 342/83 = 0.095 


By formula (16а), 


Effect of Wing on Derivative Il. 


For positive angular velocity about the Z-axis (r> 0), the 
velocity on the left wing increases and on the right decreases, 
Hence, the drag on the left wing will be greater than on the right 
wing and there appears a momerit about the Z-axis acting against : 
the rotation (fig. 14). 


The damping moment of n, from the wing is due to the induced 
and profile drags. 


Curves for the determination of n, of the wing due to the 
induced. drag are given in figure 15. In order to determine n, 
from figure 15, it ів necessary to multiply the value of пу/ ty” 
ру „2 where a, is the angle of attack of the wing computed 
from the line of zero lift in degrees, The part n. of -the wing, 


which is obtained as the result of the profile drag, is small and 
need not be taken into account’, 





PThe curves for CODI UE this i dis of Es derivative пу are 
given in reference. 5. * 
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Example, Find n, of the wing: т = 1.94, À = 6,54, а, = 99, 
In figure 15, п/а, = -0.000147, n, = -0.000147 х 92 = -0.0119. 


The total value ої the derivative n, for the airplane may be _ 
obtained from the formula | .. 


Dp = Depp Пр, |. (ав) 


Example: Find п, of the airplane a, 


99, 


The characteristics of the wing are: тп = 1.94 am A = 6.34, 


The characteristics of the tail surface аге Str/S = 0.093, in 
1%/Ъ = 0.389, Ay = 1.14, 


П. = Depp + Ny = -0.0577 - 0.0119 = -0.0696 


It must be borne іп mind that аб times the interference of the 
fuselage has a strong effect on n,. For example, for one model а 


value of n, of the tail and the fuselage was obtained equal to a D LEE LE 
third of the derivative n, for the tall only. Hence, а sufficiently 

reliable determination of n, can be obtained only бу experiment. 

Тһе theoretic&l equations for determining n, шау in individual cases 


give considerable error &nd should be considered only ав rough | КЕ 
approximations?, ! | 


Equations have been derived herein for the computation of the 
rotary derivatives, that is, the coefficients of the moments obtained 
during rotational motion of the airplane (yaw or roll). Тһе deter- 
mination of the coefficients of the forces and the moment of the 
aerodynamic forces in sideslip follows. These coefficients are called 
static derivatives. 





Early in 1940 at САНІ, work will be completed on the determina- - 
tion of the derivative n, for 10 to 15 models. From this investi- 2. = 


gation, the value of the derivative пу can be obtained for similar 
airplanes, 
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DETERMINATION OF STATIC DERIVATIVE ig 


For а positive sideslip angle В the left wing moves forward, 
the Plow under the left wing tending to lift it, that is, for 
B> O the rolling moment L is positive and 


L = LgB = Bégqsb > 0 
ig » 0 (fig. 16) 


Experiment shows the static derivative 18 тау Ре found by 
the formula 


91 al 
ів aivp" (le)y = 0 ТУ з + х s + 18 fus + let (17) 
Хе 0 


where 
y angle of dihedral of wing 
X  sweepback angle of the wing 


The characteristics of the individual terms in the equation 
will now be discussed, The magnitude (18); eo gives 1; for the 
| Х=0 
wing іп the absence ої sweepback and dihedral, The derivative 
gly - 0 arises from the change in distribution of the ciroulation 
X 50 
over the wing span as the result of the sideslip. Test curves show 
that this derivative increases with increasing angle of attack, 


This magnitude cannot be computed on account of the absence 
of test data. The plan form of the wing considerably affects the 
value of this coefficient. A rectangular wing without rounded tips 
gives a value much greater than a wing with rounded tips (reference 7). 


The effect of the dihedral of the wing is taken care of by the 


o2 | 
term У З The angle ої wing dihedral у is measured between the 


plane of the chords and the plane perpendicular to the plane of 
symmetry of the airplane and passing through the root chord (Pig. 17). 
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Experiment shows that 1; ів proportional to the dihedral М 
with a coefficient of proportionality 915/ ду, which шау be deter- 
mined from figure 18, which gives the values of the derivate 15 
for a wing without a straight midwing section. 


Example: Find 1, of the wing: т = 1.94, A= 6.54, Уз 4.259, 


8 
According to figure 18, 16/99 = 0.0133. 


ді, 
їв = У Sy = 4.25Х0.0155 = 0.0565 


For a wing with а midwing section it is necessary to subtract from 
the value obtained from figure 18 the amount 


($2) m = 0.02 Бы] (18) 


where b/b is the ratio of the span of the midwing section to the 


span ої the wing. The angle у in this case is to be measured as 
Shown іп figure 17(b), 


Example: т = 1.94, A= 6.34, у = 4.259, by /b = 0.226 


lg = 4.25 (0/0155 й 0.02(0.226)2| = 0.0523. 


Tests show that the tip of the wing greatly affects the value 
of the derivative ёв. The wing tip іп vertical projection пау ре 


such as shown in figure 19, 


From specially conducted tests (reference 7), it was found 
that if the wing tip is of the shape shown in figure 19(c), the 
effective dihedral is increased by 1 to 1.59. 


Effect of Sweepback _ 

The angle ої sweepback is measured between the lines passing 
through the focal line at 0:25 chord from the leading edge ої the 
wing and the plane perpendicular to the axis of the fuselage. The 
sweepback angle Х will be considered positive for the direction 
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Shown in figure 20. Рог a positive sweepback at а positive sideslip 
angle, the flow about the left wing ls improved and the lift force 

on the left wing is increased. As a result, an increase in the roll- 
ing moment is obtained for positive sweepback = 


X >o, lg > 0 с 


The effect ої sweepback is taken into account in equation (17) by 


і | | = 
the term Хоу. At the present time sufficient data, by use of 


which the effect of the sweepback for different wings could be 
determined, are unavailable. F, Weick (reference 8) recommends 
for rectilinear wings of aspect ratio 6 the equation 


91 
a = 0.0045 Cy (19) 


where X is the sweepback angle in degrees. 
Example. X = 2.59, Cy = 0.75(а = 99) 


From equation (19) 
1 


9? 
= = ° . 45Х0.75 =. е 82 
lg =X S 2,5(0.00 ) = 0.00822 


Effect of Position of Wing Relative to Fuselage 


The previously descrited relations give 14 for isolated wings. 


Tests show, however, a very great effect of the position of the wing 
with respect to the fuselage on the derivative ів. It was found, 


for example, that a center wing gives a value for tg approximately 


agreeing with the value for the isolated wing. A high wing increases 
the effective dihedral by 1° to 1.59 and a low wing decreases the 
effective dihedral by 2° to 5°. This change in the effective dihed- 
ге) occurs, however, at small and medium angles of attack, At large 
angles ої attack different results may be obtained, which аб the 
present time, due to the absence of sufficient test data, cannot be 
even approximately predicted, 

On &ccount of the resistance cf the flow, all projections at 
the lower surface of the wing increase ёв and conversely all pro- 


jections on the upper surface decrease lg. Іп а fixed system of 
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axes, the X-axis passes below the center of pressure of the vertical 


tail surface; therefore, the vertical tail surface will also give a 


rolling moment іп sideslip. If the angle between the line connecting 


the center of gravity with the center of pressure of the vertical 
tail surface and the velocity is denoted by a, - x 


1 бар 
18% = Ky ар % Wg” sin(x - ca) (20) 


Example, К, є 0.8, 1,/5 = 0.389, 5.2/3 = 0.095, a, = 2.05, x= 120, 


% = 99 
Ву equation 20 


bat = 0.8Х0.589Х0.093Х2.05Х0.0525 = 0.00597 


it follows from the foregoing discussion that the derivative 
18 may be determined with sufficient accuracy only by experiment. 
These equations for determining the effect of the dihedral and | 
sweepback of the wing may be used for estimating the changes in the 
airplane parameters after tunnel tests. For example, if from com- 
putation of the test data it is found that it is necessary to 
increase 16 by 0.026 and the parameters ОҒ the wing are я = 3, 


Re 7 (without mid wing) the required increase in the dihedral 
Ay is obtained from the relation Ayx0.0133 = 0.026 (the value 
0.0133 is obtained from fig. 18) whence Ay = 29. If in another 
case it is necessary (for example, from considerations of longi- 
tudinal stability) to decrease the sweepback by X = 5°, the 
change in 1, oan be found. (The computation is conducted for 


the initial state of flight of the given airplane), Let Сү = 0.6 


9? 
> = 0.00450, = 0.0045Х0.6 = 0.0027 


Alg = 0.0027Х5 = 0.0135 


Hence, decreasing the dihedral by 59, ig decreases by Alg = 0.0155. 


In order to compensate for the decrease in 18 the dihedral of the 
wing must be increased 


AyXO.0133 = 0.0155 


a 
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whence 
Ay = 19 


that is, the dihedral must be increased by 19. 


DETERMINATION OF STATIC DERIVATIVE ng 


The static derivative ng = ОС,/9В characterizes the change 


in the moment of yaw with change in the sideslip angle and depends 

mainly on the area and the shape of the vertical tail surface and 

the fuselage. The center of pressure of the fuselage for the usual 

arrangement is located ahead of the center of gravity and for this 

reason the yawing moment due to the fuselage for В > 0 ia posi- | - 
tive, that is, the fuselage is unstable. Ав may be seen from fig- 

ure 21, ng < 0. | | | 


The magnitude of the coefficient ng is determined analyti- " 
cally by the equation š — - 


| p nad 
ng = баг + Nat = K S b 7 K+ "E b S | (21) 
where | 
Se lateral &re& of fuselage 
le length of fuselage 


The coefficient Ко шау be found from figure 22 (бакеп from eo шы р yet 
reference 9) 


where 

le length of fuselage (positive) 

lag distance of center of gravity from nose of airplane (positive) 
h maximum height of fuselage (positive) 


The coefficient Кв(Кв > O) may be determined if the ratios 
1/6 ani 1 ,/1,. аге known, 
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Example: lp/h = 6.0, 1 р/Їр = 0.244 


. From figure 22 it can be seen that Кв = 0.106. 


If 82/8 = 0.55 and 16/5 = 0.595 then | 
пар = 0.106 X 0.35 X 0.595 = 0.0221. Equation (26) ів suit- 


able Гог airplanes with unstable fuselages. Sometimes the fuse- 
lage is stable (as is rarely the case) and this equation cannot 
then be applied. The value of a, may be obtained from figure 11. 


The mean value for К. is 0.8. 
Example: 2,/b = 0.589, 8.2/8 = 0.095, № = 1.14. 
From figure 11, а, = 2.05 
Пір = 7 0.8Х0.589Х0.095Х2.05 = - 0.0595 


Ng = Пар + Пар = 0.0221 - 0.0595 = - 0.0572 | ай 


Let the area of the tail surface be increased such that 
AS, p/S = 0.04. Then for the increase in the coefficient Ng 


Ang = - 0.8Х0.589Х0.04Х2.05 = - 0.0268 


Thus, when it is required to increase the vertical tail area, the 
Span must te increased but not the chord. If the vertical tail . 
surface is increased by increasing the chord, A, decreases and 


therefore a, also decreases; as а result for such an increase іп” 


бір the weathercock stability (па) із found бо undergo only 
slight increases or does not change at all. 


The equations for determining the value of ne must be сопві- 


dered to be only approximate, For а more accurate determination of 
ng, recourse must te made to experiment. As in the case of deter- 


mining the derivative ip, these approximate equations may be used 
for computing the derivative Ng after the tail area is increased 


or decreased, 
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Effect of Propeller 


In the present investigation, the effect of the propeller slip- 
Stream on the. dynamic-stability derivatives is not considered. The 
derivatives ng and n, evidently will te greatly affected by the 


propeller slipstream. Up to the present time , however, sufficiently 
accurate test data, by which а correction for the effect of the pro- 
peller on the tall surface may be introduced, are unavailable. 


DETERMINATION OF STATIC DERIVATIVE Jg 


The static derivative yg (flg. 23) is always positive 
(Jg > 0). The derivative Ур characterizes the increase іп the 
lateral force (у) acting on the airplane during an increase іп 


the sideslip angle (В). Evidently, the value of Ув depends | PER 


essentially on the size and shape of the fuselage and of the ver- 
tical tail surface. Diehl (reference 10) recommends the following 
empirical equation for computing the derivative Ув: 


ble | а | u | M mon 
Ув = 0.12 = -- | = (22) 
Example: b = 14,53, ір = 8.67, `S = 33.4 


ув = 0.12 A. 0.12 1$:95 X9.87 . | 0.453 


55.4 = 


For airplanes with usual fuselages, the mean value of уд is 0.4; 
therefore, the derivative y, has little effect cn the character- 
istics of dynamic stability. 


5. EXPERIMENTAL DETERMINATION OF DYNAMIC-STABILITY | 
DERIVATIVES 


The dynamio-stability derivatives are divided into two grcurs; 
namely, the derivatives of the moments and forces with respect tc 
the angular velocity, which are called the rotary derivatives, and 
the derivatives of the moments and forces with respect to the side- 
slip angle, which are the static derivatives, 
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The rotary derivatives шау be determined іп а wind tunnel on 
special apparatus, The procedure of these tests is, however, rather 
complicated and they will not be considered herein (reference 11). 


The static derivatives ng, ір, and ув may be readily 


obtained from the results of the usual tests in the wind tunnel. EN - 
For this purpose it is necessary to have the coefficients cf узе. 
yaw and roll and the lateral force (Cy, Су, and Су) for vari- и 


ous angles В at the given angle of attack. Curves are given in 
figures 24 vo 25 for the change іп the coefficients Cp, Cz, and Cy 


with sideslip angle. | 


ТР these coefficients are given in the fixed system of axes | 
assumed аб САНТ (fig. 1), then in determining the derivative n5, | а 
the tangent to the curve С, = f(B) for В = О must be drawn "EM 


and the slope of the angle of inclination measured (taking account 
of the scale). Special attention must be paid to the determination 
of the slope of the curve C, = Р(В). The slope of this curve i ыы 


characterizes the weathercock stability, For stable airplanes, the | 5754 
curve is like that shown in figure 24, that ів, the coefficient С, 


decreases with increase in the sideslip angle В. For neutral or | 
laterally unstable airplanes, the curves differ to a marked degree = ` pum 
from those in figure 24 near В = О only, in the range from В = -39 
to В = 459, where for neutral airplanes these curves run parallel 
to the axis of abscissas and for unstable airplanes the slope of the 
curve changes sign, that is, the coefficient С, in this range 


increases with increase in В. For а sideslip angle less than 
В = -3° or greater than В = 459, the curves in all cases are ав 
shown in figures 24, that is, the coefficient С, decreases witn 


increase in the sideslip angle В. Hence, in determining the NEA қ 
weathercock stability, the points lying within бле range В = 350 | 
must be considered. During the tests, points must be taken every 
19 angle near В = 0 but beyond |8] = 3° points may be taken 
every 29 to 50, The computations required for cbtaining the deriv- | 
ative are presented in figures 24 to 26. The coefficient 57.3 ів ы = дш 
introduced into the computation for converting degrees into radians | 
іп measuring the angle В. For а sufficiently accurate determina- ` 
tion cf the derivatives ng, lg, and ўд, it ів necessary to com- 


pute them for 4° or 5° angles of attack, (See appendix 1.) 
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Determination of Moments of Inertia of Airplans 


In order to compute the lateral dynamic stability of the air- 
plane, the coefficients of the moments of inertia, which can be 
cbtained from the following formulas, must be known; 


А АА. 424 «С 
= “То ап m 
А" wp С = qp 


where 
А polar moment of inertia with respect to X-axis 


С polar moment of inertia with respect to Z-axis 


b wing span 


m mass of airplane 

Because for various conditions of flight of the airplane (dir- 
ferent angles of attack a) the axes X and Z rotate with respect 
to the airplane, it is evident that the moments of inertia will vary 
likewise. If the momenta of inertia of the airplane about the prin- 
cipal axes (А' and С!) are known, the moments of inertia about the 
Х- and 2-вхев may be computed. by the following equations: 


А-А! сова Q 4 C' sin^9, C=C' cos? Фуд! sin? Ф, Е = А502 sin 29 


(23) 


Computations show that up to angles of attack а = 159 to 209 it 
may be assured with a sufficient degree of accuracy that A = А! 

and C = С', that ів, the moments of inertia about the X- and 
Zeaxes are equal to the corresponding moments of inertia about the 
principal axes. As will te shown, the centrifugal moment of inertia 
has little effect on the stability characterigtics and therefore Е 
need not be computed, 


At the present time no sufficiently accurate statistical for- 
mulas exist by which the moments of inertia of an airplane may be 
quickly determined and it is therefore necessary to resort to 
analytical computation in.each case. Formulas and the computing 
procedure of the moments of inertia are given in reference 12, 
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The values of i, and i, change within the ranges 
С.049 < 1, < 0.089 and 0.09< Іс“ 0.154. 


Determination of “Airplane Density" H 


The “airplane density" и is determined ty formula (11) 


where 
G/S load/sq meter, (kg/m ) 
7 unit weight of air 


b span, meters 


The computation of the dynamic stability is usually conducted for 
the altitude Н = 3000 meters and therefcre y = 0.908 kilograms 
per cubio meter, For the altitude Н =.3000 meters for standard 


atmosphere 


Example: 


Ug 


« 78 kilograms per square meter 


р = 14.55 meters 


н із а nendimensional magnitude, 


33 
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Determination об T 


The magnitude T is determined by equation (11) 


HE. 
= OSV 
In determining и for H = 5000 meters 


G 


E : ; 


where G/S is in kilograms per square meter, V is in meters per 
second, and т 18 in seconds. 


Example: G/S = 78 kilograms per square meter and У = 47.4 meters 
per second m 


78 I 
= ed 17.4 - 1.79 sec. 





4. SOLUTION OF SYSTEM CF EQUATIONS OF MCTION 


А solution of the system of equations (13) will now te sought 
in the form 


В = ве; р = pert and т = ве 


By substituting іп equations (15) and dividing by ө\, the following 
equations are obtained: 


C C 
1 L, L, | 
(ла 2 вх) зы Ран (Etone -)) В = О 
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The system of linear homogeneous equations has a solution if the 
determinant of the system is equal to zero, 


С С Е 
2.1 Др 1 
№ + 5 yp А и > > «С ten 6-1) 
1 7 і 1 
T edv £ ТА in À + in = 0 
n T 
Fe Пр ТЕ 5. СЕ 
AC зы” 7 +10 


In expanding the determinant, an equation of the fourth degree 
in A is obtained, which is called the characteristic equation. ` 


A4A* + АЗАЗ + Aga? + AV A+ AQ О ` (26) 


I | 


56 


Analysis of the coefficients snows that up to angles of attack from 
15° to 20° the magnitude i, has little effect on the roots of the 


characteristic equation and without large error ір = 0. 


If in addition, this analysis is restricted to the considera- 
tion of level flight, that is, tan @ = 0, the following equation ів 


obtained : 


А4 = 1 .. 

1 NE: 
Bx e ЛЕТ y wu 
TE = NS С Р о с A: 
часа Ца жү E (27) 

Cy 1 
Aib pol pe Ао ое 
1742 Ву ia lc іс ty in iy in ig 2 iy 
4 Ст, ftir ng ір пу 
05% 21144145 1416 
Because 


by < 0, ty 7 0, п,< 0, пр < 0, Ср? 0 


u > 0, 1470, ig > 0 


Тһе following relation is usually true: 
Аз 7 0, Ao 7 0, A4 ^ О 


The cosfficient Ay, the free term of the characteristic equation, 
may be either positive or negative. | 
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Example: 
= ya = 0.24; +В - 1.175; Be _ 0.312 
2 В T ід Й іс Й 
2 2 
T 
I=- Е go eoll 


0.37; u = 11,8 


In order to find the coefficients of the characteristic ишып, 


use is made of equations (27): 


2 . 
Az -+ УВ - v = = 0,24 - [ - 7.27 -~ 0.76] = 8.27 


2 2 І 
А To 7 do ід 2 “В in іс іс 


= {7.27X0.76 + 0.52X3.11} - 0.24 Г- 8.03] + 11.8х0.312 


1 1 
^ = (2 ) PE Pepe BE ната 


y 
Š ід io fo fa іс ta ТА ip 


= 0.24 {7.145} + 11.8X0.312X7.27 + 11.8X1.175X0.52 + 


11.8Х0.57Х1.175 = 40.809 
А КЕСЕРІ 5.11Х0.512 + 
QE Lu. ОР Ааа, Hostes 


1.175Х0.76) = - 0.5562 


= 12.7526 
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If the roots of equation (26) are denoted by Лу, Ма, ^х, and 
Ад» the solution of the system (13) is obtained in.the form ' 


ла Not Ab NC 
B = Ве * + Boe ^ + Bae ә + Bye ° 


+ Рое + Р;е + Pa e 


ню | 
T 
ду, 
r2 
Ф 


_ М АБ А А. 
Y = Re + Ree + Rze + Rye і 


Evidently, in order that the sideslip angle and the angular 
velocities of the airplane В, р, ani T should decrease with time 
(that is, іп order that the airplane should be dynamically stable), 
it is necessary that all real roots be negative and that the com- 
plex roots have negative real parts. Іп the general casé, the solu- 
tion of an algebraic equation of the fourth degree is laborious. in 
the given case, however, by making use of the special characteristics 
of the equation of lateral stability, this equation may be quickly 
and simply solved. The determination of the roots of equation (26) 
takes no more than 20 minutes, 


For airplanes of the usual type, the characteristic equation 
of the lateral stability has two real roots of which опе is very 
large and the other very small. апі two complex conjugate roots., The 
small real root is denoted by ^, the large réal root by Ag, and 


the complex conjugate roots by Az and Ад. Because ^ ів very 


small, the following equation may be written with & large degree of 
accuracy: 


A,À + Аў = O їй 
from which Aa is immediately determined!  . .. = 


eT ee зми 


lOpractical. computation shows that the value A, = 2Ag/ A4 is 
accurate, For example, if А у ів graphically determined, the 
obtained value always agrees with A, = -Ag/A,. 
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A 
Му „= - —9 28 
к (28) 


Now it can be seen that if Ag <0, then A, > О and the airplane 
is unstable. | B 


Exampie: Ау = 40.809 | Ад = -0.3362 


From formula (28) 


Ao -0.3362 
№ = - А. = ` 40.809 = 0.00825 


In determining Às, it must be recalled that À з is large and there- 
fore if the powers of А. below the third are neglected, the approxi- 
mate value is obtained 


Ма” + Ade = 0 
whence 


Ao en - Аз 


Computations on many airplanes show that the best approximation for 
Ag is the value | 


3 
А 2 ^ M . (29) 


Because і,» О and аб below-stalling angles of attack ёр - 0 
Ар < O. | 


Ву setting 


Аб) «Лё AGAS + А.Х + Ay № + Ag 


llUsually the value of ly іл gives а very small error in the 


determination of the root ^s. However, а more accurate determination 


is necessary to compute the complex roots of the characteristic equa- 
tion (more accurately, the real part). 
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Рог а more accurate determination of the root А, А ів plotted 
ав а function of A near а value equal to Л S lia, by laying 


off А оп the abscissa and А ор the ordinate. The intersection 
of the curve A(A) with the abscissa gives the accurate value of 
№. Because ` 


A'(X) = 4/3. + ЗАЛ + 2A2À + A1 
and as а polynomial for large numbers has the sign of its greatest 


term and А. is negative, usually А"(А2) < О. Hence, if after 


the first trial the point falls above the X-axis, the next value 
of A, in smaller absolute value is taken, whereas if the point 


falls below the X-axis the next value in greater absolute value is 
taken, Usually, three points are sufficient for drawing the curve 
A(A) and determining Ao. 


Example: Find No from the equation 


№ + 8,27 ЛЗ + 12.75 А2 + 40.809 А - 0.3362 = б 
Because 1% = -7.26, the value A,' = -7.0 is used for the first 
value of №2! ° 
The order of computation is clear from table І, 
The second, third, and fourth powers of A, which are required 
for the computation, may be found in appendix 2, which gives the 
powers of numbers from 2 to 15 in steps of 0.2 an3 covers the entire 


range of values of À encountered in computation, The curve A(A) 
from which Ло = -7.5 is taken, is given in figure 27. After having 


found A, and Ao, the conjugate complex roots Az, Ад are obtained 
from the following considerations: 


Equation (26) may be put in the form 
A* + AAS + ААВ + д + Ао = (A- А.а - Хе 2 +a) + b) «о 
or . 
№ + [а - (N +a JAS + - 804 +^2) + МАЮ + 
[2А - D(A, А.) A+ МА. = O 
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p By equating the coefficients of the same powers of À the following 
equation 1в obtained: 


m О 
райо А.) + МА 
у ала = С елш) 
Ао = р Ма 


whence b = A./À №; but by (28) Ау = -Ag/A,, hence | “2 ыа, 


> > 
о (N 
И ! 


p * 
її 


b E | (30) 


The value of а is found from 
Ao р - а, (АХ. Аз) + МА 
р + ^^? - Ag 


i cte (51) 


а = 


If а and b are known, the roots Az ani Ay are obtained from 
the equation 


A? ьал + b = 0 TEES 


= e - | . 
Ns a = - 2 £ 1 ъ- B P 


If № ={ + М, Ал = Ё - in anal? n | шо 





““таазе formulas are written on the assumption that b> а? [45 обо ойыс 
if b< а2/4, the formulas become ae bass 


2 
u a, а. 
Ад = у 


and in this case the airplane is unstable. If А: < 0, then Ъ<0 
and instability results. A change in the sign of А. is generally 


due to the loss of weathercock stability of the airplane, 
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5--2 (52) 


| 2 | 
1 = b - => | (55) 


Example: Find the complex roots of the equation 


AS + 8.2719 + 12.752 + 40.809 À- 0.3362 = 0 


It has been established that 


№ = + 0.00825 and Ар = - 7.5 


40.809 


- 7 ъ) = 5.61 and by equntions (31) 


A 
From equation (30), b = - E = 


БАМа - Ag _ 5.61 + 0,00825(-7,3) - 12.75 


а = XT + 0.00825 - 7.5 жы 


Tt сап be seen with the aid of equations (52) and (33) that 


t= sS). - 9-899 сазе 


nz ь- (2) = 2.315 


5. MOTION CHARACTERIZED BY ROOTS A4, Ao; А5, AND А 4 


The disturbed motion of an airplane is characterized by the 
expressions 


À a$ 5 t t 
1 № АЗ M. 
В = ве + Be + Be + В,6 
hat хаб t t 
p = Р,е + Рое + Ре + P e 
t % zë Vi 


HI 


А 4 
= Еі9 1 + вае 2 + Re + Ве 
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This motion мау be broken down іпбо three components 


" ME _ X: й 
I Ві = Ве ; P4 = Be 3 ту = Re | g Я ба : JEN 


hat A t = Nat мі = Ast A t 
ІІІ B, = Ве + Ве  , Pa = Pse +P ae , r; = Ro + Rye ` u E 


The computation of the coefficients B4, B2, Bs, By, Рі, Ро, and во 


forth will not be considered іп detail. 1% should be noted that of 7 
all these magnitudes only four are arbitrary апі are determined from 

the initial conditions, whereas the remainder can be obtained from 

these four magnitudes. 


From а numerical analysis of these ccefficients the following 
conclusions may be made: 


І, Consider the group 


+ CHE RT E 
Ві = Bet mE о 
ль _ 
== 1 
P4 = Fe ри 
_ Ж. - w 
et 22 


The magnitude Py is small in comparison with B4 and Ry. The 
root м is very small, Evidently, if the airplane is stable, 


N < O and the magnitudes. В (sideslip), p (angular velocity of 


roll), and r (angular velocity of yaw) decrease with time (with | 
increase іп Ж = tf). 
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The root À, characterizes spiral stability of the airplane, А 


Іп order to clarify the physical significance of spiral stability, 

the motion after a roll of the airplane is considered. Let the 

airplane roll in such a manner that the right wing drops. Whereas ыы Be 
in the horizontal position of the wing the. lift force balanced the 

force of gravity, the force of gravity will now give a projection 

on the Y-axis and the airplane begins to sideslip (8) on the 

right (lowered) wing. As a result of the sideslip, two aerodynamic 

moments appear: the moment 7,XxB tending to rotate the wing in 


the horizontal . position and the moment ngXP decreasing the side- 
Slip angle В. Because of the &ction of the moment Ng XB the air- 


Plane begins to yaw, that is, an angular velocity r appears, which 
gives rise to two new moments: the moment n Xr opposing the moment 


пахв and а moment 2,xr, which opposes the moment 2а ХВ. 
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Тһив, the characteristics їв and n, act favorably in elim- 
inating the disturbance and the characteristics ng and by act 
unfavorably (with respect to spiral stability). 


As can be seen from equation (28) spiral stability is charac- 
terized by Ас, which is determined by the equation 


и — 
Ag = г (2-04 ай 7,28) 


where 
M >0 520 
>0; ng «0; в.“ 0; ів > 0 


Hence, there will be stability Гог the case where Пе іп absolute 
value is greater than lyng whence і 








18| [im 
ng| |> 
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The degree of stability is characterized by the time in which the 
initial disturbance decreases by one-half 


A,t 1 
есі = 5 
— 1 
À$ = 1, 2 = - 0.695 
but 
r= d тим 
T MUS | m | 
враг = 7 0.695 у- - (34) 


In the case ої instability ^+ > 0 and the time of decrease of the 


initial disturbance by one-half will be negative. In the case of 
stability a, «0 and eer will be positive 


= 1.85 вес À. = + 0.00825 


Sapir =- 0.693 5.00825 = - 154 sec 


It has been found thet 2.5 minutes after the initial instant the 
initial disturbance will increase to twice the value. Such motion 
evidently will not бе observed by the pilot. 


Experience shows that spiral instability is not felt by the 
pilot if the time of increase of the initial disturbance to twice 
the value is large. It is to be supposed that if 

вра! > 40 - 50 вес 


at the initial flight conditions the airplane will be ss 
with regards to spiral stability, 


ТТ. Consider the group 


t 
_ Not 
ро = Рое 2 
_ Not 
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Numerical analysis shows that during a disturbance Bo and Ro are 
small in comparison with Po; that is, this group is found to char- 


acterize a motion which differs little from a pure roll. It has 
been shown that 


NS м “р 
ід 


and because at below-stalling angles ly < 0, Р will decrease very 


quickly with time because Ао is large. Generally, at the below- 
stalling angles this motion is not noticeable after 1 second, 


III, Consider the group | = 


E $ _ У К 
Pz = Fe + P,e 

- Mt АБ 

T = Rz € + R, e E 


Because А, and Nas B, and B, are conjugate complex quantities 


Аз = Р + in, Вх = bs + ib, = B(cos В + 1 sin в) х Ве iB 


A4 = $- in, By = b,- 1b = Blcos В - і віп В) _ Ben iP 


t » - in)t 
_ ве 18 „(5+ in)t oie ip АС n) 


: ttr (16 +. B): „е^ (19+ B)1] 
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But e(nt+B)1, e- (1% +В)1 -2 cos(qt- В); and, therefore, if 2B = В 


Bs = В.е + Bye = Ве cos(nt + В) 
‘Similarly, 
_ hat À $ tt D - 
р = Ре + P,e = Pe cos(nt+ y) 
Е hot hat _ 6 = 
rz = Rze + Rye = Re coa(nt + 8) 


This group characterizes the oscillatory damping of the motion. 


The time required.for decreasing the amplitude by one-half can 
now te found, | s 


tto 1 
e 22 = 
_ _— te 
Eto = - 0.693 ы dec 
Therefore 
т 
to = - 0.695 £ (35) 
Example: Т = 1.85 E = - 0.495 
| 1.83 | 
| %2 = = 0,693 -0.495 = 2,56 вес 
The period of oscillation is determined from the relations 
тт = 2% 7. 
ы» ЗД 
Шет (36) 


Example: Т = 1.83 sec п = 2.31 
1,85 85 _ 


With regard бо the oscillatory stability of the airplane, the 
following condition must be imposed: The airplane should be stable 
in oscillation and the time of decrease of the initial disturbance 
by one-half (t5) should be less than 1+ 1.5 times the oscillation 


period 


ба < 1231.5Т 


The stability in oscillation is affected mainly by the lateral static 
stability Nye The other coefficients also affect the oscillation of 


the airplane but to a smaller extent. Whereas an airplane may have a 
small static instability, it may nevertheless be dynamically stable 
but the period of oscillation (T) and the time of decrease of the 
initial disturbance by one-half (tə) will be large, At greater 


static instability, the airplane will be unstable in oscillation; 
that is, the amplitude of oscillation of such an airplane will 
increase with time. For a coefficient of lateral static stability 
of the order ng = 0.01 and greater, the least disturbance will 


teni to increase the roll of the airplane (a kind of slow spin). 
Whereas from tests in the wind tunnel of the model with locked con- 
trols it was found that the airplane is neutral or even has a small 
weathercock stability, in flight such an airplane will be unstable 
because of the mobility of the control surface, For this reason, 
from the requirement of oscillatory stability it follows that the 
airplane should be statically stable over the entire range of flight 
angles of attack and especially at below-stalling angles of attack, 


6, ORDER OF COMPUTATION OF DYNAMIC STABILITY 
Computation of Airplane Northrop 2Е for Dynamic Stability 


In this section, the entire dynamic-stability computation is 
presented from setting up tables of the structural parameters to 
obtaining the final dynamic stability characteristics. At the 
start of the computation, a table must be set up in which all para- 
meters of the airplane required for computing the lateral dynamic 
stability are included, 
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G, weight of airplane, kilograms 
Wing 

1. b, wing span, meters 

2. Drys midwing-section span, meters 

5. S, wing area, square meters 

4. Cy; chord at root, meters 

с, chord at tip, meters 

6. W, dihedral angle - angle between plane of chords and. 
plane perpendicular to plane of symmetry of airplane 

` and passing through the chord аъ the tip, degrees = 
(fig. 15) 

7. XX, sweptback angle - angle between lines passing through 
focal line taken at one-quarter chord from leading edge 
ої wing and plane perpendicular to axis of fuselage 
(X> 0 for forward sweepback), degrees (fig. 19) 

8. У, angle of twist of wing, degrees 

Vertical tail surface 


1. Sp, vertical teil surface area (area of tail with part 
of fuselage, (fig. 11), square meters 


2. Ту, tail span (measurement shown in fig, 11), meters 


5. 2+, distance from center of gravity to hinge of rudder, 
meters 


4. Angle between line of zero lift and line joining center 
| of gravity of airplane with geometric center of verti- 
cal tail surface (with part of fuselage) 
Fuselage 


Is Se, area, of projection of fuselage on plane of symmetry of 


airplane, square meters 
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50 
2, lp, length of fuselage, meters 
5. h, maximum height of projection of fuselage on plans of 
symmetry, meters 
4. Тог, distance of center of. gravity from nose of fuselage 
meters 
The magnitudes computed under I, II, III, and IV are always 
positive. | 
The table of structural parameters Рог the Northrop airplane is: 
С = 2600 kilograms Быр = 5,11 square meters 
5 = 55.4 square meters b. = 1.88 meters 
b = 14.53 meters Lh. - 5.65 meters 
bay = 3.29 meters Әр «11.69 square meters 
Co = 2.9 meters ір = 8,67 meters 
ск = 1.5 meters h = 1.45 meters 
у = 7° 50 | Ing = 2.11 meters 
X = 2° 30 x = 129 
y = 1° 50! 


The airplane is a low-wing single-engine type 


G 8 | 
g = 78 ka/n^ АЁ _ 0.095 
8 
© 1 
0 5 
À = = = 6.54 № = 1.14 
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Determination of Aerodynamic Coefficients 


l. The curve of variation of the coefficient of lift (Ст) 


against the angle of attack (a) obtained from tests in the wind 
tunnel Гог a model of the Northrop airplane is given in figure 2, 
The angle of attack is measured from the wing chord, Because Сү, 


must be determined for angles of attack measured from the line of 
zero-lift, table II is obtained. 


2. The rotary derivatives lp, lp, По, and n, were determined 


from the aforementioned equations and curves and also experimentally 
оп а special apparatus (reference 11). 


(a) Determination of 1, 


The derivative io is obtained from the data of the table of 


structural parameters (fig. 6). The curves experimentally obtained 
and by figure 6 are given in figure 50. Theoretically, the deriva- 
tive by does not depend on the angle of attack but during a test 


а certain change in ty with change in a, is usually obtained, 


The agreement obtained, as shown in figure 30, is sufficiently good 
and for computing the dynamic stability the value of by may be 


taken from figure 6 without resorting to experiment. 
- (ь) Determination of пр 


The theoretical value of np js found from figure 8 by using 


the given т and A in the table of structural parameters, The 
value of n /da taken from figure 8 must be multiplied by a, 
(angle of attack in deg. from the zero lift line), The values of 
по obtained from test ani from figure 8 are given in figure Ss 


It is obvious that а rather large disagreement results. А syste- 
matic numerical analysis of the dynamic stability equations shows, 
however, that the dynamic stability characteristics depend little 
on the velue of Dy and for this reason only an approximate value 
of np need be known, Hence, for computing the dynamic stability, 


the value, which is obtained from figure 8, may be used with suffi- 
cient accuracy. 


(c) Determination of 1, 


51 
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The determination of the derivative 1, must be divided into 
the following three stages: 


(1) Determination of 12, for an untwisted wing 


(2) Determination of the correction on 1, due to the twist 
of the wing | 


(5) Determination ої the part of 1, due to the vertical tail 
surface 


1. The derivative I, ів obtained from figure 10 by the given 
п and А of the wing. The value 1,/с, taken from figure 10 must 
be multiplied by a, (from the line of zero lift), The curve of 
1, against a, for the Northrop wing:is presented in figure 52. 
The dashed curve gives. 1, for the wing without twist. 


2, Ав the Northrop wing is &erodynamically twisted such that 
the angle of attack at the tip of the wing ів greater than the 
angle of &ttack at the tip of the midwing by 1.59, figure 12 must 
be used to determine the correction. This correction evidently 
does not depend on the angle of attack and therefore the entire 
straight line undergoes parallel displacement, The correction for 
twist is 


lp, twist = 0.0144 


5. The correction on i. due to the vertical tail. surface is 
obtained from equation (14). The required value of a, is deter- 


mined from figure 11. Тһе continuous line іп figure 52 shows the 
derivative 2. of the Northrop airplane with account taken of the 


twist of the wing and of the vertical tail surface. The figure 
also gives the points obtained from experiment. 1% is apparent that 
the agreement ів good and that the derivative l, may be computed 


theoretically without recourse to experiment, 
(à) Determination of n, 


The derivative n, ів made up of а part due to the fuselage 


with vertical tail surface and а part due to the wing. The part 
n, due to the fuselage with tail is computed from equation (16a). 
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The value a, is obtained from figure 11. The part п. due to the 


wing is computed from the data presented in figure 15. Іп order to 
obtain np, the value n/a from figure 15 must be multiplied by | 


the square of the angle of attack. The angle of attack must te in 
degrees and measured from the zero-lift line. The вит of these two 
parts gives the derivative n, of the airplane. The values of tke 


derivative n, for the Northrop airplane obtained from test and 
from computation are shown.in figure 33, I 


15 should be remembered that the derivative n, plays а great 


part in computing lateral dynamic stability and. Рог this reason it 
is desirable that n, be computed &s &ccurately &s possibie. 


Remembering also that equation (16а) may in certain cases give я 
large error, this magnitude must be experimentally determined. 


(e) Static derivatives Ур» ір» ng 


The static derivatives are determined from standard wind-tunnel 


tests; therefore these magnitudes as a rule must be obtained as shown 
in section 5. The values of yg obtained by equation (22) (curve) 


and by experiment (points) are presented in figure 34. Equation (22) 
generally gives a sufficiently accurate value of Ув. The values 

of 18 for the Northrop airplane obtained from test and from compu- 
tation are given in figure 55, The upper dashed curve gives be и 


for the value W = 7.259 (dihedral angle). Іп order to take account 
geometrically of the effect of the fuselage for a low wing in the 
formula, the dihedral must be decreased, as has been pointed out, by 
2° to 56. An airplane is considered the chassis of which is not 
streamlined and 811 projecting parts on the lower surface increase 


the effective dihedral; therefore, for the computation, the V angle | 


minus 3° is used, that is, an angle of 4,259, The lower dashed curve 
in Pigure 35 gives the value of ёв obtained from computation | 


(Че = 4.259). The value 919/9V is obtained for the given т end 
À from figure 18. | 

The effect of the midwing on the derivative is taken into 
account by the decrease in 919/9V by &n amount that may be obtained 
by use of equation (18). The center line іп figure 35 gives the 


value of їв with account taken of the part due to the vertical tail | 


surface ani the continuous line gives the value ір, if the sweep- 
back of the wing is considered, 
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The values ої ng = ас /98. obtained from test and by equa- 


tion (21) are given in figure 56, Та equation (21), the value 
К. = 0.8 is used; the value of the coefficient Ка ів obtained 


from figure 22 and a, from figure 11. 


Notwithstanding the sufficiently good agreement that was 
obtained for the Northrop airplane by computation and test for the 
values 12, = oc, /ав and ng = ӘС, /ӘВ, these formulas must be con- 


sidered only approximate and in all cases the values la and ng 
determined from the test curves C, = f(B) and С, = f- (В). 


Moments of Inertia of the Northrop Airplane 


The moments of inertia, as bas been previously stated, may be 
obtained by computation. For the Northrop airplane, the moments 
of inertia are determined by full-scale experiments according to 
the method explained in the report by U. A. Pobyedonostsev (ref- 
erence 13), The results of the investigation are as follows: 


Чл = 0.0573, іс = 0.0959 
where the values are given ln nondimensional form, 


For convenience of computation, a table of all coefficients 
that enter the characteristic equation must be set up as shown in 
table 3, in which for the derivatives the test values were taken 
and the angle a, was computed from the line of zero-lift, 


From the data of table 3, the coefficients of the character- 
istio equation Az, А, А), and A, are computed by equations (27). 
The values of these coefficients for the Northrop airplane computed 
for the angles of attack a, = 1°, 5°, 9°, and 13° are given in 
table 4 and the values A4, №, 8, and т, computed by the egua- 
tions of section 4, are given in figures 37 and 38. A comparison 
of the computed results with flight tests of reference 14 is made 
in figure 59, | 


Complete agreement was obtained for the period of oscillation T. 
The disagreement that exists for Ў (the time of decrease of the 


amplitude to 1/2) may be explained by the fact that no account мав 
taken in the computation of the effect of the propeller slipstream 


Y 
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on the vertical tail surface. The magnitudes A, and А. in 


flight have not been determined; thus it is impossible to check 
the accuracy of the computations of these values. It can be seen 


from figure 59 that the computation gives sufficiently good results, 


7. APPROXIMATE FORMULAS 


In this section, equations are given for estimating roots 
quickly from the given aerodynamio coefficients. These formulas 
may also be used in computing the corrections on Л], №, b, and m 
АҒ after computing the lateral dynamic stability it is decided to 
change the area of the vertical tail surface or the dihedral of 
the wing. 


- 


For the root Ag, the equation Аҙ = 1/44 сап be used, The 


approximate expression for A,, which determines the spiral stabil- 


ity, is obtained from Л: = -Ag/A, and because 


therefore 


Cy, 
2 (1gn, - long) 
^4 = <= — енарннаноктннно 


(37) 
| C 
ів (2 ic- ap) apiy 


The time of decrease of the initial disturbance to 1/2 is obtained _ 


by the equation 


T 


For the magnitude & (giving the time of decrease in the ampli- 
tude of oscillation tə to 1/2) the following equation can be 


used: 


(58) 
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It can be seen from equation (58) that the value cf Ё increases 
in absolute value with increase in magnitude of the damping of the 
yaw пл and the value yg, which characterizes the lateral force 


in sideslip, The magnitude of Ё. drops in absolute value with 
increase in іс ОҒ the nondimensional radius ОҒ inertia about the 


Z-axis., Vibrational instability will occur for Š > О and, as oan 
be seen from equation (38), vibrational instability may be obtained 
for а large spiral stability if А“ 0, which is usually obtained 


for large values of the derivative їв and. in the absence of 
weathercock stability, that is, when n, 20. For this reason, 


when the airplane has small weathercock stability ог instability, 
the wings should not be given a large transverse dihedral because 
the derivative ip will then be large and this condition will lead 


to, excessive spiral stability and, what is particularly important, 
may give rise to oscillatory instability, 


If the area of the vertical tail surface is changed, “the change 
in n, ів obtained by the equation 


a "S 
1.6 bá Sete i 
An, = 7 soat b S ' 


а = (9C, / да), (obtained from fig, 11) 


1+ distance from hinge of rudder to center of gravity | 


45; change in vertical tail surface area (AS> 0 for increase 
in area; AS < О Гог а decrease) 


The change in Ё is obtained from the equation А b xn, /21C, 


The new value of . Ё is obtained from Ё! = Š + АЁ, From the newly 
obtained value L', t, is found from equation (35). 


to = = 0,695 E | 5 


where T = 1.1 (G/S)/V (for an altitude of 3000 п). 
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Because A, may be approximately determined by the equation | 


Р vu 2 (а-а) аа] 
на а СҰР? 46d. 
and 


1 
№ а че n ab, where те E Ur 


iA 


the following equation is cbtained: 


i | 
n zu Ë T (2. 2) - E| (39) 


whence the period T is determined by the equation 


T = 6.28 т? where T «1.1 


Тһе greater т is the smaller the period of oscillation is 
found to be and therefore the period decreases with decrease in the 
absolute value of the derivative ng, with increase in ig, the non- 
dimensional radius of inertia with respect to the Z-axis, and with 
decrease in tne value of the derivative ёв. If the expression under 


the square root sign becomes negative (as may occur in the case where 
ng > 0, that is, in the absence of weathercock stability) the magni- 


tude т does not characterize the period of oscillation but together | 
with the magnitude € determines the asymptotic deviation of the 


airplane from its ccurse, 


The previously derived approximate equations give very good 
agreement with exact computation for the Northrop 2-Е airplane. In ` 
other cases, generally speaking, the agreement may be less “Complete 
but a qualitatively correct result will be obtainedls, 





151% the time the present article was being prepared for pub- 
lication the previously derived fcrmulas were checked Гог a large ` 
number of airplanes, Іп almost all cases the approximate formulas 
give very good agreement with the exact formulas and they may there- 


fore be used in computing the lateral dynamic stability of the EFE: 


airplans, 
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APPENDIX 1 


Up to the present time investigations in the wind tunnel have 
been conducted іп а syatem of axes different from that assumed in 
this paper ani for this reason equations are presented herein with 
the aid of which the old coefficients may be reduced to the соеї- 
ficients used in this paper. : 


The sideslip angle В (previously 7) in the diagrams of the 
wind-tunnel tests has the same sign ав В assumed in the report, 
Тһе moment coefficients have the opposite sign. Moreoyer, up to 
this time the moment coefficients were referred to рҮ in order 
to reduce them to nondimensional coefficients. If the yawing 
moment (in the old notation Сау) is based on the distance between 


the center? of gravity and the rudder hinge and the rolling moment 
(Сох) is based on the span, 


Су = ~ ёбцх Cn = - 2 р бшу Су є - аб» 


If the results from the tunnel investigations (С), Сау; C, ) 
are given in a gravity system of coordinates 9С,,/9В, p 
and 90,/98 аге obtained from the curves and the coefficients lg, 
ns and yg are obtained from the equations . 


д0, OC x : 
bp 798 = - 2 sg Х57.5 


ac ль, ӘС әс, ӘС 
= 38” = - а Х57.5; Ур = 38 =- 2 wg Х57.5 


If the results are given іп a system of chord axes, all moments 
and forces are based on ove, the rolling moment C',, is based. on 


the span, and the moment of. yaw С uy is based on ty (the distance 
between the center of gravity and the hinge of the rudder), 
+ | | ty 


бау = С”, 008 9 + Ь Су біп аз Cy = Су 208 a- Cu sin а 


PT LE - 


2074. 
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and finally 


OC, OO rs | OC "ух 1% дошу 
lg 7 dB 272795 57.5 = 7 el 35 cos a + B "55 sina 57.5 


APPENDIX 2 
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SECOND, THIRD AND FOURTH POWERS OF NUMBERS FROM 2.0 TO 15.0 


CÓ со DO го ro 
" е е е ә 
(со Ф н [Ó O 


Q Q CI Q G: 
Со б WS 


8.0 
8.2 
8.4 
8.6 
8,8 


я» = Н № 
Q ол нь Оо 


Ф 


сл <л ол J (л 
р о d> NOS 


219,0|3242,0|47978,6 


15.0/225.0/3375.0 [50625.0 


о о о 
A го O 





ұзат Wh VOVN 


ы | 
) 2074 аа ' 
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TABLE ТТТ - COEFFICIENTS FOR COMPUTING 


DYNAMIC STABILITY OF NORTHROP AIRPLANE 


ЕСЕЈЕ аач Е 


Ув 0855 0.55 0.48 0.45 
ТВ 0.058 0.088 0,068 0.068 
ng -0.053 -0.057 -0.03 -0,032 
Ly -0.45 - 0,45 -0.42 -0.58 
ік | 0.042 | 0.115 | 0.180 0.26 
ny -0.019 -0.054 -0.05 -0.065 
п, |-0.055 -| -0.0675 |-0.075 | -0.089 

1 1 - 

a. DLL PEE 4 

in s iles’ ic - 10.4 

TABLE IIIA 
| |as |а = 50 |а = зе |o = 15° | 
Ст, 
y 
= 0.265 0.265 0.24 0.225 
y: 
B 

YA ШОН TS T: 00267 1.175 ЕТӘ 
n 
i -0.396 | -0.385 | -0.312 -0.555 
1 
T -7.3 -7.3 «Т 27 -Б.57 
"А 
ly | 
чу 0.727 1.99 5,11 4.5 
= -0.197 -0.354 - , 052 -0.676 
Пу 
іс -0.572 -0.702 -0.76 -0.717 


Н = 11,8 


64 


12.75 | 40.809 
15.32 | 44,14 
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Figure 5. - Projection of gravity force on axis during roll of airplane. 
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Figure 6, - Value of derivative 
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Гог wings ої various aspect ratios and tapers. 
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Figure 8. - Values of n/a, for wings of various aspect ratios and tapers, . _ 
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Figure 10. - Values of 72,/са for wings.of various aspect ratios and tapers. 
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Figure 12. - Graphs for determining correction on 2, due to twist of wing. 
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Figure 13. 
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- Determination of correction on 77, for twist of wing, 
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Figure 15, - Values ої п/а for wings of various aspect ratios and tapers. | 


та МАСА ТМ 1264 





Figure 16. - Derivative Lae 
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Figure 17, - Measurement of dihedral angle of wing (deg). 
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Figure 18. - Curves Гог determining derivative 18 ОҒ wing with dihedral, 
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21. = Derivative рв. 
Figure š 
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- Curves for ый effect ої fuselage on derivative пд. 


Figure 22, 





- Derivative Ув. 


Figure 25. 
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Figure 24, - Variation of увчіпе moment coefficient С, with sideslip angle В, 
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Figure 25. - Variation of rolling moment coetficient G 
"m OC, (Сз — 10° 57.3 0,0125.57.3 


2 with sideslip angle В. 


{ — — = — ШШЕ a — n €—M——M € 
8 38 i0 0.071 


76 


МАСА ТМ 1264 





© де 3 
== f =05 





Figure 27, - Determination of root Ao mo MT equation. 
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Figure 09, - Curves of lift force and drag against angle of attack. 
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Figure 30. - Derivative ір ої Northrop airplane. 
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Figure 51. - Derivative п, of Northrop airplane. 
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Figure 52. - Derivative 7, of Northrop airplane. 
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Figure 55, - Derivative п, of Northrop airplane. 
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Figure 56. - Derivative Dg = дс. /9в. 
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Figure 57. - Values of roots of characteristic equation X. and Хо (Northrop airplane). 


81 


82 NACA ТМ 1264 






DX 







ву | | |_ НЕЕ Еж 
ТЕ Lis 












TP ЕШ И ШШ ЖШ ИЕ 
LLL 
Ae В А ЕЕ ыы 
| N |l ||| Ë ° a ба» from flight twat 





eta | | | T=— 


EHE Rana n s sus — 
XE вё Эй 
HE Li [T+ | 


oe bs : Е 


Figure 39. - Comparison of values T and t, obtained by computation Sind from flight test 
(Northrop airplane}. 
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